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Abstract 

This  report  contains  a compilation  of  molecular  dynamics  computer  simulation  results 
for  the  simple  point  charge  extended  model  (SPC/E)  for  water.  In  addition  to  published 
results,  this  compilation  contains  unpublished  results  developed  by  the  authors.  The  prop- 
erties included  are  the  static  dielectric  constant,  the  self- diffusion  coefficient,  the  internal 
energy,  and  the  pressure.  The  simulation  values  of  the  dielectric  constant  are  compared 
with  a new  correlation  based  on  available  experimental  values. 

Introduction 

The  simple  point  charge  extended  (SPC/E)  model1  for  water  has  proven  to  be  one  of  the 
more  robust,  but  simple  of  the  many  interaction  models  proposed  for  water  since  the  early 
work  of  Stillinger  and  Rahman.2  In  this  note  we  present  a survey  of  SPC/E  results  for 
the  potential  energy  U,  the  self-diffusion  coefficient  Ds,  and  the  dielectric  constant  e over 
a wide  range  of  temperatures  and  densities.  No  claim  for  completeness  of  this  survey  is 
made,  nor  for  the  validity  of  any  specific  result.  There  is  reasonably  good  consistency 
between  the  results  in  various  studies.  This  report  is  the  outgrowth  of  some  simulations 
performed  by  the  authors.  Also,  experimental  values  of  the  dielectric  constant  are  listed 
so  that  the  extent  to  which  this  model  mimics  water  can  be  examined. 

The  SPC/E  model  assumes  that  a water  molecule  is  a rigid  object  with  an  intramolecular 
distance  of  0.1  nm  between  oxygen  and  hydrogen  interaction  sites  (O-H  distance)  and  with 
an  angle  of  109.47°  between  the  O-H  bonds.  The  intermolecular  site-site  interactions  are 
defined  in  terms  of  the  distances  between  the  sites.  Although  these  sites  are  commonly 
interpreted  in  terms  of  oxygen  and  hydrogen  atoms,  they  are  properly  only  sites  for  atom- 
atom  and  Coulomb  interactions.  There  are  partial  charges  assigned  to  the  sites  so  as  to 
mimic  an  effective  charge  distribution  of  a water  molecule  in  liquid  water.  The  charge 
on  the  oxygen  site  is  -0.8476e  and  the  charge  on  the  hydrogen  site  is  0.4238e  with  the 
result  that  the  permanent  dipole  moment  fi  — 7.84  x 10~30C  m (2.35D).  The  resulting 
charge  distribution  has  a quadrupole  moment  (in  units  of  10~4OC  m2)  QXI  of  7.290, 
of  -6.286  and  Q:z  of  -1.007  compared  with  the  experimental  values  of  8.77,  -8.34,  and  -0.43 
respectively.3 
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The  permanent  dipole  moment  is  larger  than  the  moment  of  an  isolated  molecule,  /i o = 
6.17  x 10-30C  m (1.85D),  and  the  induction  energy  associated  with  the  excess  dipole 
moment, 

Eind  — /^o) 

j 

where  a is  the  polarizability  of  the  molecule,  amounts  to  5.213  kJ/mol.  Inclusion  of  this 
mean-field  value  of  the  induction  energy  when  the  site  charges  were  determined  is  the 
distinction  between  the  original  SPC  model  and  the  SPC/E  model. 

In  addition  to  the  Couloumb  interaction  between  molecules,  the  model  has  a Lennard-Jones 
interaction, 

<f>(r)  = 4 eLJ[(a/r)12  - (a /r)6} 

between  the  oxygen  sites.  The  Lennard-Jones  parameters  are  eu  = 78.24 kg  and  a — 
0.31656  nm.  Here  Icb  is  Boltzmann’s  constant.  The  energy  of  the  system  is  the  sum  of  the 
mean-field  induction  energy,  the  Coulomb  energy,  and  the  Lennard-Jones  energy. 

First,  the  unpublished  simulation  results  of  the  authors  will  be  presented.  Then,  results 
of  other  studies  will  be  listed.  A summary  discussion  concludes  this  report. 

Molecular  dynamics  studies  of  SPC/E  water 

A series  of  molecular  dynamics  simulations  have  been  performed  by  RDM  using  the  SPC/E 
model  to  determine  the  dielectric  constant  of  this  model  at  elevated  temperatures  and 
at  densities  below  1000  kg/m3.  The  development  of  an  extensive  correlation4  of  ex- 
perimentally determined  dielectric  constants  in  this  range  of  thermodynamic  states  was 
the  motivating  factor  for  this  study.  This  correlation  is  based  on  data  in  the  interval 
230 K < T < 873 1\  with  densities  on  the  order  of  1000  kg/m3  or  less.  This  correlation 
reproduces  the  experimental  values  within  ±0.5  with  much  narrower  bounds  for  temper- 
atures below  the  critical  temperature.  Our  comparisons  between  simulation  values  and 
values  provided  by  this  correlation  are  limited  to  this  range  of  temperatures  and  densities. 
If  no  correlation  value  is  quoted,  it  is  for  this  reason.  If  a value  for  the  dielectric  constant 
is  obtained  by  an  extrapolation  outside  the  range  of  the  data,  that  value  is  enclosed  in 
parentheses,  (...). 

The  simulations  were  performed  on  systems  of  216  rigid  water  molecules  contained  in 
a cubic,  periodically  continued  cell.  Ewald  summations5  with  conducting  conditions  at 
infinity  were  utilized  to  account  for  the  long  range  Coulomb  interactions  while  the  Lennard- 
Jones  interactions  were  truncated  at  one-half  the  box  size.  The  equations  of  motion  were 
integrated  using  an  iterated  Beeman  algorithm6  and  a time  step  of  10~15  s.  A quaternion 
representation'  was  used  for  the  orientational  degrees  of  freedom.  The  temperature  of 
the  system  was  controlled  by  separate  Nose-Hoover  thermostats8  for  the  translational  and 
rotational  degrees  of  freedom.  The  coupling  constant  between  the  thermostats  and  these 
degrees  of  freedom  had  a value  Q — 0.05mcr2  where  m is  the  mass  of  a water  molecule  and 
a is  the  Lennard-Jones  parameter.  This  value  provided  a reasonably  steady  value  for  both 
the  translational  and  rotational  kinetic  energies  as  the  simulation  progressed. 

The  properties  reported  here  are  the  energy,  self- diffusion  coefficient,  and  the  dielectric 
constant.  A total  of  sixteen  states  have  been  simulated.  The  temperature  and  density 
pairs  along  with  these  properties  are  listed  in  Table  1. 
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The  duration  of  a given  simulation  was  dictated  by  the  long  sampling  times  needed  to 
obtain  stable  estimates  for  the  dielectric  constant,  e,  from  the  fluctuation  expression5, 

4tt  < M2  > 

+ T kBTV  ' 

Here,  < M2  > is  the  ensemble  average  of  the  square  of  the  total  dipole  moment  of  the 
system  of  molecules,  V is  the  volume,  and  T is  the  absolute  temperature.  Since  each 
configuration  yields  only  one  sample  for  M and  the  decay  of  the  correlations  of  the  total 
moment  is  a relatively  slow  process,  requiring  tens  of  picoseconds,  with  the  result  that  the 
simulation  times  needed  range  from  100  ps-1000  ps.9  The  time  step  used  in  the  integration 
of  the  equations  of  motion  is  10-15  s.  Thus  several  hundred  thousand  iterations  of  the 
equations  of  motion  are  necessary  for  each  thermodynamic  state. 

The  values  of  the  dielectric  constants  listed  in  Table  1 have  converged  to  within  ±1  unit. 
This  conclusion  is  based  on  the  fluctuations  in  the  long  averaging  time  estimates  of  < M2  > 
with  increasing  averaging  time.  The  absence  of  any  apparent  drift  in  the  average  was  taken 
as  a necessary  condition  for  convergence.  The  self-diffusion  coefficients  were  estimated 
using  the  Einstein  relation  for  the  long-time  behavior  of  the  mean-square  displacement  of 
molecules, 

< [Ar]2  >— > 6Dst 

where 

< [Ar]?  > = < i - r j ( 0 )] 2 > 

3 = 1 

and  the  angle  brackets  < . . . > indicate  an  average  over  the  enclosed  quantity. 

A second  set  of  simulations  was  performed  in  the  supercritical  region  by  AW.  These  simula- 
tions employed  216  rigid  water  molecules  enclosed  in  a cubic,  periodic  cell  under  constant 
volume  conditions.  The  long-range  electrostatic  interactions  were  include  via  an  Ewald 
simulation,  whereas  the  Lennard-Jones  interactions  were  spherically  truncated  at  0.93  nm. 
The  equation  of  motions  were  integrated  using  a velocity  Verlet  algorithm  coupled  to  a 
SHAKE/RATTLE  procedure0  to  maintain  the  intramolecular  water  geometry.  A timestep 
of  2.0  fs  was  used,  velocities  were  periodically  rescaled  to  maintain  a constant  temperature 
and  the  translational  and  rotational  temperatures  were  monitored  to  avoid  any  imbalances 
between  these  degrees  of  freedom.  After  a 20  ps  equilibration  period  the  systems  were 
sampled  for  an  additional  400  ps  during  which  the  ensemble  averages  were  calculated.  The 
results  from  these  simulations  are  listed  in  Table  2.  No  uncertainty  estimates  have  been 
made  for  these  results. 

Results  from  other  simulation  studies  of  SPC/E  water 

The  properties  of  the  SPC/E  model  have  been  examined  for  a wide  range  of  densities  and 
temperatures.  A systematic  examination  of  the  density  and  temperature  variation  of  the 
dielectric  constant  of  water  has  been  reported  recently.10  The  results  of  this  study  are 
listed  in  Table  3.  The  results  of  other  studies  for  the  dielectric  constant,  for  the  potential 
energy  (including  the  induction  energy  Eind ),  and  for  the  self- diffusion  coefficient  are  listed 
in  Tables  4 and  5,  along  with  the  references  for  the  source  of  the  values.  Table  4 contains 
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dielectric  constant  results,  Table  5 contains  potential  energy  results.  Diffusion  coefficient 
results  are  found  in  Tables  4 and  5. 

The  results  from  Ref.  13  lie  along  the  liquid  vapor  coexistence  curve  for  the  SPC/E  model. 
The  temperature-density  trajectory  determined  there  does  not  exactly  coincide  with  that 
of  water,  with  the  largest  deviation  of  the  model  density  on  the  liquid  side  being  about 
15%  below  that  of  the  liquid  near  the  critical  point.  The  deviations  on  the  vapor  side  are 
large  with  the  model  density  greater  than  the  density  of  water  vapor. 

The  long-range  Coulomb  interaction  has  been  treated  in  various  ways.  In  the  simulations 
leading  to  Tables  1 and  2,  Ewald  summation  methods  were  used.  This  is  generally  agreed  to 
be  a “correct”  method,  although  it  is  computationally  expensive.  Alternative  formulations, 
such  as  reaction-field  methods  can  yield  good  results  if  used  judiciously.  Truncation,  with 
or  without  smoothing,  tends  not  to  produce  the  proper  structure  in  the  liquid.  This  is 
particularly  troublesome  for  dielectric  constant  determinations. 

In  addition  to  results  summarized  in  Tables  3,  4,  and  5,  the  surface  tension  of  SPC/E  water 
has  been  determined  with  reasonably  good  agreement  with  experimental  values  between 
316K  and  573K.19.  There  is  also  one  report  of  the  shear  viscosity  at  277K  and  at  300K.16 

Discussion 

Since  the  simulation  results  listed  in  Table  1 were  motivated  by  interest  in  the  high  tem- 
perature values  for  the  dielectric  constant,  it  is  appropriate  to  examine  that  body  of  results 
in  more  detail.  The  density-temperature  values  for  which  SPC/E  values  of  e have  been 
reported  are  indicated  in  Fig.  1.  The  temperature  variation  of  the  dielectric  constant  is 
indicated  in  Fig.  2 where  a dashed  line  joins  values  on  the  coexistence  curve  reported 
in  Ref.  13  and  solid  lines  join  points  of  constant  density.  A logarithmic  scale  is  used 
to  avoid  visual  clutter.  There  is  fairly  good  consistency  between  the  values  obtained  by 
various  workers  for  states  with  temperatures  greater  than  400K.  There  are  however,  some 
considerable  discrepancies  reported  for  lower  temperatures  (and  necessarily  higher  densi- 
ties). The  discrepancies  occur  in  the  region  where  one  must  be  careful  to  obtain  adequate 
sampling  for  the  total  moment  fluctuations,  (M2).  In  general,  the  SPC/E  values  for  the 
dielectric  constant  are  lower  than  the  experimental  values  obtained  at  the  same  density 
and  temperature.  This  is  indicated  in  Fig.  3 where  the  relative  deviation  of  the  simulated 
value  of  the  dielectric  constant  from  the  value  provided  by  the  data  correlation  of  Ref.  4. 
Figs.  4 and  5 display  the  relative  deviations  along  isochores  and  isotherms,  respectively. 

The  self-diffusion  coefficients  and  potential  energy  values  for  the  SPC/E  model  are  collected 
in  “one  place”.  No  attempt  has  been  made  to  compare  the  calculated  and  experimental 
values  for  these  quantities. 

All  of  the  results  reported  here  are  based  on  classical  dynamics  simulations.  Since  the 
de  Broglie  wavelength  of  a hydrogen  atom  at  300K  is  on  the  order  of  0.1  nm,  quantum 
dynamical  effects  can  be  expected  to  be  important.  This  has  been  shown  to  be  the  case 
for  room  temperature.20  It  is  noteworthy  that  the  isotopic  mass  effect  on  the  self- diffusion 
coefficient  for  hydrogen,  deuterium,  and  tritium12  is  smaller  than  the  quantum  dynamical 
effect. 
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T,  K 


Fig.l  The  density  and  temper- 
ature values  for  which  the  di- 
electric-constant of  the  SPC/E 
model  has  been  reported.  The 
symbols  indicate  the  source: 
square,  Table  1;  circle,  Table  2; 
*,  Ref.  13;  open  triangle,  Ref. 
14;  solid  triangle,  Ref  10;,  and 
diamond,  Refs.  11  and  12. 


T , K 


Fig  2.  The  temperature  de- 
pendence of  the  SPC/E  di- 
electric constant  is  shown. 
The  symbols  have  the  same 
meaning  as  in  Fig.  1.  A 
dashed  line  joins  the  points 
on  the  SPC/E  coexistence 
curve.13  Solid  lines  join  con- 
stant density  points  from  a 
single  source.  Starting  from 
the  bottom  of  the  plot,  the 
solid  lines  mark  densities  of 
200,  300,  400,  600,  800,  and 
1000kg/m3. 
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Ae/e  Ae/e 


Fig  3.  The  relative  deviations 
of  the  simulation  values  from 
the  correlated  experimental 
values  of  the  dielectric  constant 
are  displayed.  The  symbols 
have  the  same  meaning  as  in 
Fig.  1. 


Fig  4.  The  relative  deviations 
of  the  simulation  values  from 
the  correlated  experimental 
values  of  the  dielectric  constant 
are  displayed  for  isochores,  as 
noted  in  the  legend.  The  sym- 
bols have  the  same  meaning  as 
in  Fig.  1. 
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Fig  5.  The  relative  deviations 
of  the  simulation  values  from 
the  correlated  experimental 
values  of  the  dielectric  constant 
are  displayed  for  isotherms,  as 
indicated  in  the  legend.  The 
symbols  have  the  same  mean- 
ing as  in  Fig.  1. 


Table  1.  The  dielectric  constants,  the  energy,  and  the  self- diffusion  coefficients  for  SPC/E 
water  produced  by  RDM  are  listed  for  a set  of  temperatures  and  densities,  along  with  the 
values  of  the  dielectric  constant  of  water  provided  by  the  correlation.4  Values  in  (...)  are 
extrapolated  from  the  correlation. 


T,K 

Pm,  kg/m3 

e,  SPC/E 

€,  corr. 

U,  kJ/m 

Ds,  10  9 m2/s 

430 

1000. 

43.5 

48.6 

-34.8 

10.2 

573 

720. 

20.8 

20.5 

-25.8 

31.8 

573 

800. 

23.4 

24.0 

-27.4 

27.0 

573 

862. 

26.4 

26.9 

-28.4 

23.8 

573 

920. 

27.3 

29.7 

-29.4 

21.8 

573 

1000. 

31.3 

(33.7) 

-30.4 

18.1 

673 

800. 

18.8 

19.7 

-24.6 

32.7 

673 

862. 

20.4 

22.1 

-25.8 

28.8 

673 

1000. 

24.0 

(27.8) 

-27.8 

21.2 

823 

800. 

14.6 

15.7 

-21.4 

40.5 

823 

862. 

15.9 

(17.6) 

-22.5 

34.5 

823 

1000. 

19.1 

(22.3) 

-24.5 

26.9 

873 

540. 

7.4 

(8.2) 

-15.1 

65.8 

873 

720. 

12.1 

(12.5) 

-19.0 

49.2 

873 

900. 

15.9 

(17.7) 

-22.2 

34.8 

878 

858. 

13.7 

(16.3) 

-21.5 

38.0 
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Table  2.  The  dielectric  constants,  the  energy,  and  the  pressure  for  SPC/E  water  produced 
by  AW  are  listed  with  the  appropriate  temperatures  and  densities.  Also  listed  are  the 
values  of  the  dielectric  constant  of  water  provided  by  the  correlation.4 


T,K 

Pm,  kg/m3 

e,  SPC/E 

e,  corr. 

U,  kJ/m 

P,  MPa 

673 

300. 

4.7 

4.7 

-18.4 

31.5 

773 

100. 

1.8 

1.7 

- 7.6 

25.0 

773 

200. 

2.8 

2.7 

-12.4 

41.5 

773 

300. 

4.1 

4.1 

-15.8 

57.0 

773 

400. 

5.7 

6.0 

-18.6 

79.0 

773 

500. 

7.6 

8.3 

-21.1 

111.5 

773 

600. 

9.6 

10.8 

-23.4 

166.0 

773 

700. 

12.3 

13.7 

-25.4 

258.0 

773 

800. 

14.4 

16.8 

-27.5 

421.0 

773 

900. 

17.3 

-29.3 

672.0 

773 

1000. 

19.9 

-30.7 

1053.0 

873 

300. 

3.7 

3.7 

-13.9 

84.0 

973 

300. 

3.4 

-12.5 

109.0 

1073 

300. 

3.1 

-11.4 

135.0 

1173 

300. 

2.9 

-10.5 

159.0 

1273 

300. 

2.7 

- 9.7 

184.0 
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Table  3.  The  dielectric  constants,  the  energy,  and  the  pressure  for  SPC/E  water  determined 
in  Ref.  10  are  listed  for  a set  of  temperatures  and  densities.  Values  of  the  dielectric  constant 
of  water  provided  by  the  correlation  are  limited  to  500K.4 


T,K 

Pm,  kg/m3 

€,  SPC/E 

e,  corr. 

500 

800. 

28.0 

28.8 

500 

1000. 

38.6 

40.0 

1000 

200. 

2.5 

1000 

400. 

4.8 

1000 

600. 

7.7 

1000 

800. 

11.0 

1000 

1000. 

14.7 

1500 

200. 

1.9 

1500 

400. 

3.2 

1500 

600. 

4.8 

1500 

800. 

6.6 

1500 

1000. 

8.5 

2000 

200. 

1.7 

2000 

400. 

2.5 

2000 

600. 

3.5 

2000 

800. 

4.7 

2000 

1000. 

6.1 

2500 

200. 

1.5 

2500 

400. 

2.1 

2500 

600. 

2.9 

2500 

800. 

3.8 

2500 

1000. 

4.7 
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Table  4.  The  dielectric  constants  and  self-diffusion  coefficients  quoted  in  various  studies 
of  the  SPC/E  model.  No  entry  for  Ds  indicates  that  no  value  was  reported.  No  entry  for 
e,  corr  indicates  that  the  state  point  is  outside  the  range  of  the  correlation.4 


T,K 

Pm,  kg/m3 

€,  SPC/E 

e,  corr. 

Ds , 10-9  m2/s 

Reference 

263 

998. 

112.  ±16. 

92.1 

0.96 

11 

298 

998. 

89.  ±13. 

78.5 

2.15 

11 

298 

998. 

71. 0±  6.1 

78.5 

12 

373 

958. 

51.  ± 4. 

55.5 

6.8 

11 

300 

998. 

83.1 

77.9 

2.6 

13 

373 

949. 

66.1 

54.9 

7.9 

13 

473 

841. 

30.0 

33.8 

19.6 

13 

570 

36. 

1.4 

556.0 

13 

570 

674. 

18.7 

18.7 

34.7 

13 

610 

82. 

1.8 

1.7 

238.0 

13 

610 

545. 

14.3 

12.3 

53.6 

13 

620 

112. 

2.6 

2.0 

191.0 

13 

620 

507. 

12.4 

10.8 

54.1 

13 

630 

162. 

3.2 

2.6 

157.0 

13 

630 

470. 

9.6 

9.5 

65.6 

13 

640 

217. 

4.1 

3.4 

129.0 

13 

640 

426. 

9.4 

8.0 

71.7 

13 

652 

326. 

6.0 

5.4 

106.0 

13 

324 

1000. 

72.0±14. 

70.7 

14 

377 

1111. 

55. 7±  4.5 

14 

518 

1000. 

39. 1±  3.0 

38.2 

14 

772 

876. 

18. 7±  1.1 

19.4 

14 

782 

1000. 

21. 1±  1.1 

14 

810 

257. 

3.9±  0.2 

3.3 

14 

878 

858. 

15. 1±  0.9 

14 

975 

756. 

10.0±  0.5 

14 

994 

500. 

6.1±  0.3 

14 

994 

909. 

14. 2±  0.9 

14 

1040 

1000. 

13. 3±  0.8 

14 

1042 

429. 

5.2±  0.2 

14 

1091 

702. 

8.8±  0.5 

14 

1158 

808. 

9.3±  0.3 

14 

1278 

1000. 

10. 0±  0.6 

14 

11 


Table  5.  The  potential  energy  and  self-diffusion  coefficients  quoted  in  various  studies  of 
the  SPC/E  model.  No  entry  for  U or  for  Ds  indicates  that  no  value  was  reported. 


T,K 

Pm,  kg/m3 

U,  kJ/mol 

Ds,  10  9 m2/s 

Reference 

210 

1024. 

-48.01 

0.05 

15 

220 

1025. 

-46.95 

0.16 

15 

230 

1027. 

-46.08 

0.31 

15 

240 

1026. 

-45.08 

0.50 

15 

250 

1025. 

-44.19 

0.59 

15 

261 

1022. 

-43.20 

1.08 

15 

282 

1016. 

-41.28 

1.51 

15 

307 

1001. 

-39.60 

2.51 

15 

328 

987. 

-37.26 

3.57 

15 

263 

998. 

-43.6 

0.96 

11 

298 

998. 

-41.4 

2.15 

11 

298 

998. 

-41.5 

12 

373 

958. 

-37.1 

6.8 

11 

277 

990. 

1.75 

16 

300 

976. 

3.18 

16 

306 

998. 

-41.4 

2.5 

1 

273 

900. 

1.36 

17 

273 

1000. 

1.46 

17 

300 

800. 

2.92 

17 

300 

900. 

2.70 

17 

300 

1000. 

2.43 

17 

611 

500. 

-20.37 

18 

996 

500. 

-11.90 

18 

1979 

500. 

- 4.65 

18 

294 

1000. 

-41.96 

18 

1038 

1000. 

-21.26 

18 

1986 

1000. 

- 9.86 

18 

951 

1500. 

-21.68 

18 
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